Layer-by-layer entangled spin-orbital texture of the topological surface state in Bi2Se3 
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With their spin-helical metallic surface state, topological insulators (TIs) define a 
new class of materials with strong application potential in spintronics. Technological 
exploitation depends on the degree of spin polarization of the topological surface state 
(TSS), assumed to be 100% in phenomenological models. Yet in Bi2Se3, an archetypical 
TI material, spin- and angle-resolved photoemission spectroscopy (S-ARPES) detects a 
spin polarization ranging from 20 to 85%, a striking discrepancy which undermines the 
applicability of real TIs. Here we show —studying Bi2Se3 by polarization-dependent 
ARPES and density-functional theory slab calculations — that the TSS Dirac fermions 
are characterized by a layer-dependent entangled spin-orbital texture, which becomes 
apparent through quantum interference effects. This explicitly solves the puzzle of 
the TSS spin polarization in S-ARPES, and suggests how 100% spin polarization of 
photoelectrons and photocurrents can be achieved and manipulated in TI devices by 
using linearly polarized light. 



Topological insulators (TIs) define a new state of mat- 
ter in which strong spin-orbit interaction (SOI) leads 
to the emergence of a metallic topological surface state 
(TSS) formed by spin-nondegenerate Dirac fermions jl]- 
[6]. To capture the physics of TIs, a spin- momentum 
locking with 100% spin polarization is usually assumed 
for the TSS in time-reversal invariant models [3]-[5]. A 
fully spin-polarized TSS has tantalizing properties, such 
as robust edge modes and exotic quasiparticle excita- 
tions, making TIs good candidates for spintronic appli- 
cations [7]; and also a playground for engineering Majo- 
rana fermions, and exploiting their non-Abelian statistics 
in topological quantum computation [8]. The successful 
realization of topological insulating behavior in quantum 
wells m [To] and crystalline materials pTHT3] brings us 
closer to the practical implementation of theoretical con- 
cepts built upon novel topological properties. 

To this end, the fundamental question is that of the 
effective spin-polarization of the TSS Dirac fermions in 
real materials and devices, i.e. how realistic is the hy- 
pothesis of a 100% spin polarization. The large discrep- 
ancy in the degree of TSS spin polarization determined 
for Bi2Se3 by spin-resolved ARPES (angle-resolved pho- 
toemission spectroscopy) - ranging from 20 to 85% [M]- 
[18j - highlights the complexity of real TIs. First princi- 
ple density-functional theory (DFT) calculations indicate 
that the TSS spin polarization in members of the Bi2X3 
material family (X=Se, Te) can be substantially reduced 
from 100% [I9l[20]. In addition, based on general sym- 
metry arguments, it was shown that the spin polarization 
direction of photoelectrons in spin-resolved ARPES can 
be very different from that of the TSS wavefunction [21]. 
However, the role played by the intrinsic properties of 



the TSS wavefunction in defining the highest spin polar- 
ization that could be achieved, for instance in d.c. and 
photoinduced electrical currents, has remained elusive. 

We report here that the TSS many- layer- deep exten- 
sion into the material's bulk - in concert with strong SOI 
- gives rise to a layer-dependent, entangled spin-orbital 
texture of the Dirac fermions in Bi2Se3. This enables 
the precise control of both in- and out-of-plane spin po- 
larization of the photocurrent in spin-resolved ARPES - 
all the way from to ±100% - by varying energy, po- 
larization, and angle of incidence of the incoming pho- 
tons. More generally, the layer- by- layer variation in spin- 
orbital entanglement is of fundamental importance to 
spintronic studies and applications, whose outcome will 
depend critically on how one couples to the TSS. A re- 
markable consequence, which we will specifically exploit 
in this study, is that one can gain exquisite sensitivity to 
the internal structure of the TSS wavefunction, ^tss, via 
quantum interference effects using ARPES. In particular, 
the spin-orbital texture is captured directly in the linear- 
polarization dependence of the ARPES intensity maps in 
momentum space, and can be fully resolved with the aid 
of ab-initio DFT slab-calculations. 

Let us start our discussion from the Bi2Se3 ARPES re- 
sults in Fig.[l] measured with a and tt linearly-polarized 
21.2 eV photons [22]. Based on our experimental geom- 
etry (Fig.jlJ^) and the photoemission selection rules, a- 
polarized light probes mainly the in-plane Px and Py or- 
bitals, whereas 7r-polarized light a combination of both 
in-plane and out-of plane {pz) orbitals: the observed 80% 
reduction in overall intensity by switching from tt- to a- 
polarization indicates that the TSS has a dominant Pz 
orbital character. As for the evolution of the ARPES in- 
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Fig. 1. Linear polarization dependence of the measured ARPES intensity in momentum space. (A) Schematics 
of the experimental geometry, with tt (horizontal) and a (vertical) linear polarizations, and horizontal photoelectron emission 
plane. (B) ARPES dispersion measured along K— F— K with tt polarization; the zero of energy - here and throughout the paper 
- has been set at the Dirac point (DP) for convenience. (C-D) Constant energy ARPES maps from above (0.1 and 0.2 eV) 
and below (-0.1 eV, inset) the DP, measured with a polarization; (F-G) same for tt polarization. (E,H) Normalized variation 
of the a- (E) and vr-polarization (H) ARPES intensity, along the Dirac contours, plotted as a function of the in-plane angle Lp. 



tensity around the Dirac cone, in a-polarization (Fig.[T]3 
and E) we observe a twofold pattern at both 0.1 and 
0.2 eV binding energy above the Dirac point (DP), con- 
sistent with a previous report [23 , although somewhat 
asymmetric with respect to the ky = {) plane (see in partic- 
ular Fig.[l^); this suggests a tangential alignment of the 
in-plane Px,y orbit als with respect to the Dirac constant- 
energy contours. Conversely in 7r-polarization (Fig.[l^- 
H) we observe an unprecedented onefold, strongly asym- 
metric pattern at 0.1 eV above the DP, which evolves into 
a threefold pattern while still retaining some asymmetry 
at 0.2 eV; this is in stark contrast with the uniform distri- 
bution of intensity along the Dirac contour expected for 
the dominant Pz orbit als. Finally, at — O.leV below the 
DP, a threefold pattern is observed for both polarizations 
(inset in C and F). 

The asymmetry in ARPES intensity between iky is 
particularly evident in 7r-polarization at 0.1 eV in Fig.[lJ^ 
and in the band dispersion of Fig. [l^. This finding, which 
might seem in conflict with the time-reversal invariance of 
the TSS, provides fundamental clues on the structure of 
^TSS- Time-reversal invariance requires the state at +k 
with (pseudo) spin up to be degenerate with the state 
at — k with (pseudo) spin down, i.e. to have the same 
real-orbital occupation numbers. This so-called Kramers 
degeneracy, together with the ARPES selection rules for 
linearly polarized light, forbids intensity patterns which 
are different at ±k. We emphasize here that this re- 
striction can be rephrased in terms of purely in-plane 
momentum coordinates, i.e. iky, only for a perfect 2- 
dimensional TSS with a delta-function-like density, for 
which kz plays no role. Thus the observation of an imbal- 
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Fig. 2. Calculated ARPES polarization dependence. 

(A,B) Constant energy maps of the cr-polarization ARPES in- 
tensity from the topological surface state (TSS, 0.1 and 0.2 eV 
[24 ) and bulk valence band (BVB, -O.leV in the inset); (C) 
corresponding variation of the ARPES intensity as a function 
of the in-plane angle Lp. (D-F) Same data as in (A-C), but 
now for TT-polarized light. (G-I) Calculated constant-energy 
circular dichroism ARPES patterns at 0.1, 0.2, and 0.3 eV 
above the DP; the insets show the patterns obtained by ro- 
tating the sample by 90 ° about the sample normal. 
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Fig. 3. Layer-by- layer entangled spin-orbital texture of the TSS. (A) Electronic dispersion from our 250-layer-slab 
DFT model [22], with TSS in orange and bulk states in green (^tss is composed of those states that exceed a 20% projection 
onto the top 5 atomic layers in real space). (B) Percentage contribution of pz and px,y orbitals to ^tss at 0.15 eV above 
the DP, resolved layer- by- layer, for the top 15 atomic layers of the slab (Sel is the natural cleavage plane). (D-F) Layer- 
and orbital-projected charge density along the 0.15 eV /c-space contour indicated in (C): surfaces are defined by riO^cj)) — 
r' ^^(^5 0)(ttl r fc^2,T',fc)5 where i and r are the layer and local orbital basis indexes, and Z the cubic harmonics; 

the surface color represents the expectation value of the Sy operator. While the pz channel is layer independent (D), a strong 
layer-dependent spin-orbital entanglement is observed for px,y (E). The total layer-resolved TSS texture (F) is obtained by 
adding all p orbital contributions according to their relative, layer- dependent weight from panel (B). 



ance in ARPES intensity at iky , together with the estab- 
lished time-reversal invariance of TIs, necessarily implies 
that ^TSS must have a finite extent - albeit not a disper- 
sion [25] - along the third dimension. While details will 
become clear when discussing our DFT results in Fig.|3) 
we anticipate that this - together with the strong SOI 
- leads to a complex layer- dependent spin-orbital entan- 
glement in Bi2Se3, which becomes apparent in ARPES 
through photoelectron interference effects. 

Before delving into the analysis of the layer-resolved 
^TSS we note that, by performing ARPES intensity cal- 
culations [22l [26] for the extended TSS and bulk wave- 
functions from our DFT slab-calculations, we accurately 
reproduce polarization-dependent ARPES experiments. 
As shown in Fig.[2]A.-F, we obtain very different inten- 
sities at ±k|| in excellent agreement with the results 
for both (J and tt polarizations. Specifically, we repro- 
duce the approximately twofold pattern in a polariza- 
tion, stemming from the spatial configuration of Px,y oi"- 
bitals (Fig.|2]A. and B); and also the threefold pattern 
away from the DP (Fig.^E and insets of A and D), which 



originates from the hybridization between TSS and bulk 
states (and is associated with the kz dispersion of the 
latter [27]). As a final test of the robustness of our DFT 
analysis of ^tss, we have calculated constant-energy cir- 
cular dichroism ARPES patterns, which are also in ex- 
cellent agreement with previous studies |28l [29] . 

To gain a microscopic understanding of the properties 
of ^xss we present our DFT results for a 250-layer slab 
of Bi2Se3 [22] in Fig.jsjA., with bulk states in green and 
TSS in orange. The in- and out-of-plane p orbital pro- 
jections in Fig.[3p confirm that ^tss indeed has a large 
Pz (70%) character, although p^^y (30%) is also signifi- 
cant, and most importantly that ^tss extends deep into 
the solid. Even though the orbital weight decays expo- 
nentially with the distance from the surface, as expected 
for a surface bound-state, ^tss extends approximately 2 
quintupole layers (QL) below the surface (~ 2nm), with 
'^75% contribution from the 1^* QL and '^25% from the 
2^^ QL. Note also the interesting layer dependence of the 
orbital character: while for most layers the main compo- 
nent is the out-of-plane Pz^ for the 5^^ (for example) the 
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Fig. 4. TSS spin texture and photoelectron spin polarization. (A-D) Spin texture of the Dirac cone upper branch in 
Bi2Se3 obtained from the expectation value of the layer- integrated, orbital-projected spin operators, normalized to the orbital 
occupation [22]; in-plane and out-of-plane spin components are represented by arrows and colors, respectively. Note that (A-B) 
and (C-D) have different color scales but that the arrow scaling remains the same, with the largest arrows representing full 
polarization; also, moving away from F corresponds to moving along the Dirac dispersion away from the DP (^0.4eV at maps 
edge). (E-F) Prediction for the photoelectron spin polarization (P) measured in spin-resolved ARPES as a function of photon 
energy and incidence angle [22^; two experimental geometries are examined in tt polarization for the same k point located at 
0.15 eV along f — M, as indicated in the sketches [in (E) only the P component along y is shown]. 



in-plane ^x,y is actually dominant. 

As a consequence of the relativistic SOI, which directly 
connects orbital to spin flips via the /^s^ terms of the 
spin-orbit operator 1-s = /^-^^ + + /~s+)/2, the 

strongly layer-dependent orbital occupation becomes en- 
tangled with the spin polarization of ^tss- To visualize 
this entanglement, in Fig|3jl)-F we present the layer- and 
orbital-projected charge density along the 0.15 eV Dirac 
contour indicated in FigjSp, with surface colors repre- 
senting the expectation value of the Sy operator [22] . 
The p2-pi"ojected charge density, being associated with 
a single orbital, cannot be entangled and has the layer- 
independent spin helicity shown in Fig.jsjD [30]. In con- 
trast, a strong layer-dependent spin-orbital entanglement 
is observed for "p^^y because the eigenstates can be a lin- 
ear combination of Vx,^^ Py,i^ and similar states, result- 
ing in a complex set of charge-density surfaces. These 
surfaces show two overall spatial configurations having 
opposite spin helicity, which are oriented tangentially 
and radially with respect to the Dirac cone contour [31 , 
as seen in Fig.[3|E. In Fig.|3j^ we show the total layer- 
dependent charge density obtained by adding in- and out- 
of-plane contributions according to their relative weights 
in Fig. [3^; from this it is clear that while the Pz orbit als 



dominate, the in-plane Px^y orbitals lead to a substantial 
spin-orbital entanglement of the combined ^tss- 

This entanglement also leads to complex in- and out- 
of-plane spin-texture, as shown in Fig.[4]A.-D where the 
layer-integrated spin patterns of individual and total p 
orbitals are presented. While for p^ we find the in-plane 
helical spin texture expected for the TSS this is not the 
case for the Px and Py orbitals, which exhibit patterns 
opposite to one another. Combining all contributions 
[(5^*°*^V^totai in Fig.ji^], the TSS out-of-plane spin tex- 
ture vanishes in the vicinity of the DP; most important, 
the in-plane spin polarization is reduced from 100% to 
75% at the DP, and to 60% at 0.4 eV above the DP [22]. 
Note that this is also critically dependent on the relative 
Px^y orbital content of ^tss, which increases from 25% 
to 45% over the same energy range [22]. 

We will now highlight the interplay between the spin- 
orbital texture of ^tss, photoelectron interference, and 
the measured ARPES intensity. While a complete deriva- 
tion - inclusive of selection rules, Bi and Se cross sec- 
tions, as well as photoelectron escape depth - is given 
in the supporting online material [22], for the purpose 
of this discussion we approximate the ARPES intensity 
as / oc Ke^^ '^IA-pl^TSs)!^, expressed in terms of plane- 
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wave photoelectron final states for simplicity. By writ- 
ing ^TSS as linear combination of the layer- dependent 
eigenstates, ^tss = cr ^^V^i^ky with i and a being 
layer and spin indexes, the ARPES intensity then be- 
comes /ocEalEie-^''^'^(e'^i'""|A.p|a,V^-k|,)l'- Here 
the e~*^^^^ phase term accounts for the photoelectron 
optical path difference stemming from the TSS finite ex- 
tent into the bulk. Because both e~*^^^* and ip^^ vary 
from layer to layer (the latter via the relative orbital con- 
tent as shown in Fig.|3^), the photoemission intensity is 
dominated by interference between the V^^j^^^ eigenstates, 
and can in fact be regarded as the Fourier transform of 
the layer-dependent ^tss- We also note that, because 
the phase of photoelectrons is defined by additive kz and 
k|| contributions, reversing the sign of either kz or ky 
will change the ARPES intensity, i.e. I{kz) y^I{—kz) and 
especially /(ky) ky) as observed experimentally. 

Photoelectron interference also severely affects the spin 
polarization Px,y,z = ftl'l'l^lll'l^l measured in spin- 
resolved ARPES [22 . This exhibits a strong dependence 
on photon energy, polarization, and angle of incidence, 
which in general prevents the straightforward experimen- 
tal determination of the intrinsic spin-texture of Bi2Se3. 
While comprehensive results are presented in Fig. S2 and 
S3 [22 , in Fig. [4^ and F we show as an example the same 
k point along F — M measured in two different geome- 
tries, probing selectively Py^z (E) and Px,z (F) orbitals. 
In Fig.|4^, because (SPy) and (Sp^) (the spin polarization 
of the Py and Pz orbitals) are antiparallel at this specific 
k point, Py varies between ±100% upon changing 6>, and 
oscillates wildly as a function of photon energy (with the 
exception of ° and 90 ° , which probe Py and Pz sepa- 
rately). However, if the sample is rotated by 90° as in 
Fig. Up {SP"") becomes parallel to (5^^) and the measured 
Px,y,z are all independent of photon energy and incidence 
angle, allowing the detection of the intrinsic spin polar- 
ization. We note that this behavior is consistent with 
reported spin-resolved ARPES results [18]: for the sit- 
uation of Fig.|4^, Py > 80% was obtained, close to our 
100% expectation; along f — K, Py was observed to vary 
from 25% at hu = 36 eV to -50% at 70 eV, while we 
obtain +20±10% and -40±15%, respectively [22 . 

In conclusion, we have shown that the spatial extent 
of the TSS into the solid and its strong layer-dependent 
spin-orbital entanglement are responsible - via photo- 
electron interference - for the apparent time-reversal 
symmetry breaking in ARPES and the large discrep- 
ancy in the estimated TSS spin-polarization from spin- 
resolved ARPES. This is of critical importance for many 
applications and fundamental studies of TIs. For in- 
stance we note that the photoelectron interference re- 
sponsible for /(k||) 7^/(— k||) in ARPES also provides an 
explanation for the so-far puzzling observation of spin- 
polarized electrical currents photoinduced by linearly- 
polarized light ^32j, which also is associated with an im- 



balance in the number of photoelectrons removed at iky . 
In addition, exploiting photoelectron interference in spin- 
resolved ARPES provides a way not only to probe the 
intrinsic spin texture of TIs, but also - and most impor- 
tantly - to continuously manipulate the spin polarization 
of photoelectrons and photocurrents all the way from 
to ±100% by an appropriate choice of photon energy, 
polarization, and angle of incidence. 
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